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Abstract 

We present the results of measurements of two-particle angular correlations between trigger particles reconstructed 
in the ALICE Eorward Muon Spectrometer (-4.0 < ij < -2.5) and associated particles reconstructed in the cen¬ 
tral barrel detectors (|? 7 | < 1.0) in p-Pb collisions at a nucleon-nucleon center-of-mass energy of 5.02 TeV. At low 
transverse momentum (pt), the reconstructed trigger particles mainly originate from weak decays of primary pions 
and kaons, whereas at high pj, they originate mainly from the decay of heavy-flavor particles. The ridge structure 
elongated in the 77 -direction, discovered recently in p-Pb collisions, is found to persist to the pseudorapidity ranges 
studied here. The second-order Eourier coefficients of muons are extracted after subtracting the correlations obtained 
in low-multiplicity events from those in high-multiplicity events. The ratio of coefficients in the Pb-going and p- 
going directions is presented as a function of pj, and the coefficients are observed to have a similar pj dependence 
in both directions with the Pb-going coefficients larger by 16+6%. The results are compared with calculations from a 
multi-phase transport model. 
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1. Introduction 

Two-particle angular correlations in Sip and A 77 , 
where and A?; denote differences in the azimuthal 
angle {ip) and the pseudorapidity ( 77 ) between trigger and 
associated particles, have been widely used to explore 
particle production mechanisms in the field of relativis¬ 
tic heavy ion physics. In nucleus-nucleus collisions, 
ridge-like structures on the near (at Ap « 0 ) and away 
side (at Ap =» n) along the A 77 axis are generally consid¬ 
ered to be essential evidence of collective behavior of 
the created medium JT]. Recently, these ridge-like struc¬ 
tures were also observed in high-multiplicity proton- 
proton 121 and proton-lead (p-Pb) collisions 12] [H 0 
at LHC energies. These observations have opened new 
debates on their origin. Various theoretical attempts 
have been made to explain these results, such as mod¬ 
els based on hydrodynamic flow, color glass condensate, 
and multi-parton interactions. 


One way to further understand these ridge-like struc¬ 
tures is to investigate their 77 -dependence in p-Pb colli¬ 
sions. The CMS Collaboration has reported preliminary 
results evaluating the azimuthal harmonic coefficients 
from particle correlations within I 77 I < 2 , and results in¬ 
dicate a mild 77 -dependence within their acceptance 0 . 
In addition, predictions based on a 3+1 dimensional vis¬ 
cous hydrodynamical model and a multi-phase transport 
model (AMPT) at large 77 (2.5 < I 77 I < 4) show a strong 
77 -dependence in the harmonic coefficients Q . 

ALICE has measured angular correlations between 
trigger particles from the Eorward Muon Spectrometer 
(EMS) and associated particles from the central barrel 
detectors, which cover -4.0 < 77 < -2.5 and I 77 I < 1.0, 
respectively, in p-Pb collisions at yfs^ - 5.02 TeV. 
The analysis procedures and results are described in the 
following sections, while more details can be found in 
Ref. 11 . 
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2. Experimental setup and analysis 

The current analysis is based on the p-Pb collision 
data collected by the ALICE detector in 2013 with a 
4 TeV proton beam and a 1.58 TeV/nucleon lead ion 
beam. Beams were provided in both configurations, 
denoted by p-Pb (Pb-p) for proton beam going in the 
direction of negative (positive) rj, allowing to measure 
muons in the p- and Pb-going directions. Due to the 
asymmetric beam energies, the nucleon-nucleon center- 
of-mass system is shifted in rapidity by 0.465 in the 
direction of the proton beam. In order to ensure con¬ 
sistency of the multiplicity classes in both beam con¬ 
figurations, only symmetric parts of the VO detectors 
(2.8 <11 <3.9 in VO-A and -3.7 < 77 < -2.7 in VO-C) 
were used for the event classification. 

The analysis followed a similar procedure as in a pre¬ 
vious measurement 0 to evaluate v^{2PC,sub), which 
is the second-order Fourier coefficient of the muons ex¬ 
tracted from two-particle correlations with the subtrac¬ 
tion method. The associated yield per trigger particle 
was measured in a given event class using muon tracks 
as trigger particles and inclusive charged hadrons in 
the central barrel detectors as the associated particles. 
The yield from the highest multiplicity class (0-20%) 
was subtracted by the yield from the lowest multiplicity 
class (60-100%) to isolate the long range correlations, 
and fitted with 

ao+2 ai cos(A^)-i-2 02 cos(2Ai,c))-i-2 03 cos(3A^) .(1) 

The relative modulation, y„A{2PC, sub), is given by 
where b is the baseline of the low-multiplicity 
class (60-100%) estimated from the integral of the per- 
trigger yield around the minimum. v^{2PC,sub) are 
then obtained as 

v^{2PC, sub) = V„a{2PC, sub) / ,^y^"^{2PC, sub), (2) 

where y^^{ 2 PC, sub) is measured by correlations among 
inclusive charged hadrons in the central barrel detectors. 

Muon tracks used as trigger particles in the analysis 
were selected according to the same selection criteria 
as in Ref. IfTOll in the range of 0.5 < py < 4.0 GeV/c. 
The front absorber of the FMS is designed to suppress 
the contribution of weak decays from light hadrons, 
and affects the parent particle composition accordingly. 
Figure [T] shows the composition of parent particles of 
reconstructed muons in the p-going (indicated by p-Pb) 
and Pb-going (Pb-p) directions. This composition as 
a function of reconstructed py has been obtained with 
the DPMJET Monte Carlo (MC) event generator ifTTl 
and GEANT3 IfT^ for the detector response simulation. 
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Figure 1: Parent particle composition for reconstructed muon tracks 
in the p-going (p-Pb) and Pb-going (Pb-p) directions from detector 
response simulation of the ALICE Forward Muon Spectrometer with 
DPMJET Monte Carlo events. 

While weak decays of pions and kaons dominate the low 
Pt region, their contributions are significantly reduced 
at pt > 2 GeV/c and decays of heavy-flavor particles 
are dominant in this region. The parent particle compo¬ 
sitions in the p-going and Pb-going directions are simi¬ 
lar in this simulation. 

Both muon-track and muon-tracklet correlations 
were measured using central-barrel tracks and Silicon 
Pixel Detector (SPD) tracklets, repsectively, as asso¬ 
ciated particles. Central-barrel tracks, reconstructed 
based on the information from the Inner Tracking Sys¬ 
tem (ITS) and the Time Projection Chamber (TPC), 
were selected in the same way as in Ref. |[3l with the 
fiducial region \t]\ < 1 and 0.5 < pj < 4.0 GeV/c. SPD 
tracklets were reconstructed based on the primary ver¬ 
tex position and two hits on the SPD layers with It/I < 1. 
Two hits were required to have an azimuthal angle dif¬ 
ference less than 5 mrad. This requirement suppresses 
fake and secondary track contributions and caused the 
mean py of selected tracklets to be about 0.75 GeV/c, 
which was estimated using the DPMJET MC. 

While v^{2PC, sub) from muon-tracklet correlations 
was measured for both p-going and Pb-going direc¬ 
tions, v^{2PC,sub) from muon-track correlations was 
only calculated for the p-going direction due to low 
statistics. More details on ALICE detectors, such as the 
TPC, ITS, FMS and VO, can be found in Ref. 1T3]. 

3. Results 

An example of the associated yield per trigger parti¬ 
cle is shown in Fig.[^ for p-going muon tracks as trigger 
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Figure 2: An example of the (A^, A^) distribution of the associated 
yield per trigger particle in the highest multiplicity class using trigger 
particles from the FMS and associated particles from the central barrel 
detectors in p-Pb collisions. 


particles and SPD tracklets as associated particles in the 
0-20% event class. As the Arj region is far away from 
the near-side jet-peak region around (Atf, Arf) = (0,0), 
no near-side jet contribution is observed. A near-side 
ridge structure is visible in this high multiplicity class. 

v^{2PC, sub) for p-going direction was measured with 
both tracks and tracklets as associated particles, sepa¬ 
rately, and is shown in Fig. Although the two differ¬ 
ent associated particle classes probe different pj ranges, 
the v^{2PC, sub) results agree with each other through¬ 
out the whole py ranges studied here. This agreement 
can be considered as evidence of the factorization be¬ 
tween trigger and associate V 2 . 

Using tracklets as associated particles, v^{2PC, sub) 
for the p-going and Pb-going directions was measured 
and plotted in Fig. The Pb-going v^{2PC, sub) is 
observed to be larger than the p-going v^{2PC, sub) 
throughout the measured py range, and both have a sim¬ 
ilar py dependence. v^{2PC, sub) from AMPT model 
calculations are also provided in Fig. It has to be 
noted that the predictions in Ref. Q based on a 3-I-1 
dimensional viscous hydrodynamical model and the 
AMPT model cannot be directly compared with the 
data, although they cover the same rf range. Trigger 
particles used in the predictions are primary particles, 
and V 2 is extracted for primary particles while the pre¬ 
sented v^{2PC, sub) here is for decay particles after the 
front absorber of the ALICE FMS. This involves smear¬ 
ing effects in tp, and passing through the front absorber 
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Figure 3: v’^{2PC, sub| in the p-going direction, calculated for muons 
in -4 < ;; < -2.5 from two-particle correlations with the subtraction 
method, extracted from muon-track and muon-teacklet correlations in 
p-Pb. Systematic and statistical errors are plotted as boxes and verti¬ 
cal eiTor bars, respectively. 


dramatically changes the parent particle composition for 
the muons relative to the particle composition of pri¬ 
mary particles. Instead, the AMPT calculations in Fig.j^ 
were prepared for the comparison with the data. Using 
the same parameters as in ||71, generator level primary 
particles were simulated to decay into muons using the 
PYTHIA decayer lfT4l . The effects of the front absorber 
were accounted for by applying the relative efficiencies, 
which were evaluated for muons from pion and kaon de¬ 
cays relative to that for heavy flavor decay muons from 
a detector simulation of the ALICE FMS. The relative 
efficiencies can be found in Ref. iH. Other details to 
calculate v^{2PC, sub) in AMPT were the same as those 
in the data analysis, including event classihcation based 
on the number of charged particles in 2.8 < rf < 3.9 and 
-3.7 <!)< -2.7. 

The results from the AMPT calculation qualitatively 
agree with the data for py < 1.5 GeV/c, the region dom¬ 
inated by weak decays of pions and kaons as shown in 
Fig.B However, a qualitative discrepancy between data 
and AMPT calculations is observed for py > 2.0 GeV/c. 
As the contribution from heavy-flavor decays are the 
largest in this region and heavy-flavor V 2 in AMPT cal¬ 
culation is equal to 0, possible scenarios to explain this 
discrepancy can be either a hnite value for the V 2 of the 
muons from heavy-flavor decays or large differences in 
particle composition and V 2 of each species between the 
AMPT generator and data. 

To quantify the difference between v^{2PC, sub) for 
the p-going direction and v^{2PC, sub) for the Pb-going 
direction, the ratio of the two is calculated and shown in 
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Figure 4: ii^{2PC, sub| from muon-tracklet correlations in the p- 
going and Pb-going directions for -4 < j] < -2.5 in p-Pb collisions at 
a center-of-mass energy of 5.02 TeV. Systematic and statistical errors 
are plotted as boxes and vertical error bars, respectively. Results from 
AMPT calculations are compared to the data. Only systematic eiTors 
ai‘e provided for the AMPT results. 

Fig.H] Considering the size of systematic and statistical 
errors, the ratio is observed to be almost independent of 
Pj. a constant fit to the ratio is performed and returns 
the value, 1.16+0.06, with a;^f^/NDF=0.4. Results from 
AMPT calculations show a similar py dependence. 

All systematic uncertainties for Fig. [30 and are 
described in Ref. (H. 

4, Summary 

The V 2 of muons reconstructed in the ALICE Forward 
Muon Spectrometer was measured from by two-particle 
correlations with subtraction method in p-Pb collisions 
at a nucleon-nucleon center-of-mass energy of 5.02 TeV. 
With Monte Carlo simulations, the composition of par¬ 
ent particles of reconstructed muons, largely affected by 
the front absorber, was determined to vary depending 
on Pj. v^{ 2PC, sub) for the Pb-going direction was ob¬ 
served to be larger than v^{2PC, sub) in the p-going di¬ 
rection by 16 ± 6%, and both showed similar pj depen¬ 
dence. Comparison with AMPT calculations revealed a 
disagreement with the data for pj > 2.0 GeV/c, where 
muons from heavy-flavor decays dominate. 


Figure 5: The ratio of \'^{2PC, sub} in the Pb-going direction and 
\'^{2PC, sub} in the p-going direction. The ratio from AMPT calcula¬ 
tions is compared to the data. 
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